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Abstract 

r-^. \ We calculate the isosinglet unpolarized quark- and antiquark distributions at a 

low normalization point in the large-_/V c limit. The nucleon is described as a self- 
consistent soliton solution of the effective chiral theory. The ultraviolet cutoff is 
implemented by a Pauli-Villars subtraction. The quark and antiquark distributions 
satisfy the momentum as well as the baryon number sum rule, corresponding to 
"constituent" quarks whose structure is not yet resolved at the low scale. Already 
at this level a sizable number of antiquarks is present. 
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An essential ingredient in the theory of deep-inelastic scattering and other hard processes 
are the leading-twist parton distribution functions of the nucleon []1|, |2| . Their evolution in 
the asymptotic region is governed by perturbative QCD and well understood. However, 
the initial conditions for the perturbative evolution, i.e., the parton distributions at a 
relatively low normalization point, belong to the field of non-perturbative physics, and at 
present one has to rely on model calculations to estimate them. 

Recently, an approach has been formulated to calculate the leading-twist parton dis- 
tributions in the large-iV c limit, where the nucleon can be described as a chiral soliton 
@, IJ. It was shown that this field-theoretic approach allows to preserve all general re- 
quirements on parton distributions such as positivity and the partonic sum rules which 
hold in QCD. Results for the leading quark- and antiquark distribution function in the 
l/-/V c -expansion (isosinglet unpolarized and isovector polarized) show reasonable agree- 
ment with the parametrizations of the data at a low normalization point || ; the approach 
has been extended to transverse polarized distributions 0. First calculations of the sub- 
leading distributions have also been reported 0. 

The description of the nucleon as a chiral soliton is based on the dynamical breaking 
of chiral symmetry, which can be encoded in an effective theory, valid at low energies. 
The effective theory can be derived from the instanton vacuum of QCD, in which case 
the ultraviolet cutoff is determined by the inverse average instanton size, p _1 = 600 MeV, 
the scale at which the dynamically generated "constituent" quark mass drops to zero 
0. In actual calculations the fall-off of the dynamical quark mass is usually simulated 
by applying some external UV regularization. Hadronic observables, such as the nucleon 
mass, axial coupling constant etc. are, as a rule, insensitive to the details of the UV cutoff. 
When one turns to the calculation of parton distribution functions, however, there are very 
strong restrictions on how one should introduce the UV cutoff in the effective theory. The 
point is that one has to preserve certain general properties of parton distributions such as 
positivity, sum rules etc., which can easily be violated by an arbitrary UV regularization. 
Specifically, the regularization should preserve the completeness of the set of quark single- 
particle wave functions in the soliton pion field. One possible regularization method which 
fulfills all requirements is a Pauli-Villars subtraction || [|] . 

An important aspect which has not been taken into account in the calculation of parton 
distribution functions so far, related to the constraints on the UV regularization, is the 
self-consistent description of the chiral soliton. The classical pion field describing the 
nucleon should be determined as a minimum of the static energy, suitably regularized in a 
way that ensures consistency with the regularization applied to the parton distributions. 
The point where the two regularizations meet is the momentum sum rule for the isosinglet 



distribution, which is satisfied only if the classical pion field is a solution of the equations 
of motion of the effective chiral theory ||. The calculations of parton distributions in 
refs. 0, U were performed with a particular variant of the Pauli-Villars regularization and 
using a fixed soliton profile, and the momentum sum rule was satisfied only approximately. 

In the present letter we compute the isosinglet unpolarized distributions in the ap- 
proach formulated in refs.[Q, |J with a self-consistent description of the nucleon. In fact, 
a class of regularizations exist which preserve all general requirements on parton distri- 
butions (positivity, sum rules etc.) and simultaneously lead to a stable minimum of the 
energy With such a regularization the momentum sum rule is explicitly preserved (along 
with other sum rules). We also investigate the effects of possible finite regularizations on 
the valence and, particularly, the antiquark distribution. 

The calculation of parton distributions from the effective chiral theory relies on the 
parametric smallness of the ratio of the dynamical quark mass, M, to the UV cutoff, A. 
[In the instanton vacuum this ratio, (Mp) 2 , is proportional to the small packing fraction 
of instantons J7[.] The basic expressions for the twist-2 parton distributions in refs.[Q, [| 
were derived in the limit M/A —>■ 0. Differences between various UV regularizations are 
parametrically of order (M/A) 2 . Since the existence of a minimum of the energy depends 
on the details of the UV regularization, i.e., on terms of order (M/A) 2 , one does not, from 
a principal point of view, achieve a higher accuracy of the parton distributions by making 
the description of the soliton self-consistent. Nevertheless, since the existence of a stable 
soliton is an important qualitative issue - - not only for the momentum sum rule - - it 
makes sense to "improve" the regularization at level (M/A) 2 in order to ensure soliton 
stability, and this is the point of view we take here. 

Self- consistent soliton with Pauli-Villars regularization. The classical pion field char- 
acterizing the nucleon in the effective chiral theory at large N c is of "hedgehog" form 

U(x) = exp [z(n ■ r)P(r)] = cos P(r) + z(n • r) sin P{r) 

x 

r = x, n = — . (1) 

r 

The profile function, P(r), is determined by minimizing the static energy of the pion field, 

E tot [U] = N c E lev + N c Y, (En-E^), (2) 

neg. cont. 

where E n are the quark single-particle energies, given as the eigenvalues of the Dirac 
Hamiltonian in the background pion field (M is the dynamical quark mass), 

H<$> n = E n $ n , H = -i-y^dk + M-flT*, 

U^( x ) = ^^U{x) + ^—^-U\x). (3) 



The spectrum of H includes a discrete bound-state level, whose energy is denoted by E\ ev , 
as well as the positive and negative Dirac continuum, which are polarized by the presence 
of the pion field. The sum in Eq.(|2|) runs over all occupied quark single-particle levels in 
the soliton, that is, the bound-state level and the negative continuum. The EJ® denote 
the energy levels of the vacuum Hamiltonian given by Eq.(|J) with U — 1. 

The contribution of the Dirac continuum to the energy, Eq.(^), is logarithmically di- 
vergent and requires UV regularization. We consider here regularizations of the form of a 
Pauli-Villars subtraction, which can be implemented consistently with the Pauli-Villars 
regularization of the parton distribution functions below. A possible regularization of 
Eq.(H) is to subtract from the continuum contribution a suitable multiple of the corre- 
sponding sum computed with the quark mass replace by a regulator mass, M PV : 

_ 71^2 _ 

E tot , rcg [U] = N c E lcv + N c £ (£„-£(°))-iV 53 (£ PV> -4°)g,(4) 

neg. cont. PV neg. cont. 

where E P v, n are the eigenvalues of the Hamiltonian, Eq.(||), with M replaced by M P y and 
the same pion field. The subtraction coefficient follows from the fact that the logarithmic 
divergence of the unregularized sum, Eq.(fJ), is proportional to M 2 . The value of the 
regulator mass can be fixed from the (also logarithmically divergent) pion decay constant, 

F2 = w /** m» m» r d ' k M - 



(2vr) 4 (M 2 + k 2 ) 2 C M 2 PV J (2vr) 4 (M PV + k 2 ) 2 

N C M 2 M 2 PV 
= ^ l0g M^' (5) 

E n = 93 MeV. The regulator mass now plays the role of the UV cutoff of the theory, A. 
It was shown in ref. that the energy functional regularized according to Eq.@ 
possesses a minimum with respect to the profile function, i.e., that a stable soliton exists. 
The minimum is determined by the stationarity condition, 

SE t , 



7 tot,reg 

5P{r) 



sin P(r)S(r) + cos P(r)V(r) = 0, (6) 



where 



M 2 
S(r) = N c S lev (r) + N c £ S n (r) - N c — — £ <S PVj „(r), (7) 

neg. cont. 1V1 P V neg. cont. 

M 2 
V(r) = Nc V lcv (r) + N c £ V n (r) - N c — — £ V PV , n (r), (8) 

neg. cont. v PV neg. cont. 



S n (r) = MJ d 3 x$t(x)<5(|x|-r)7 $ n (x), (9) 

V n (r) = M|rf 3 x^(x)(5(|x|-r)7V(n-T)$„(x), (10) 



and similarly for Spv,n(r),Vpv,n(r). The profile function can be found by iterative so- 
lution of Eq.(|6|), see ref. ||, where also a number of hadronic nucleon observables have 
been calculated with this regularization. In the leading order of the l/iV c -expansion 
the nucleon mass is given simply by the value of the energy, Eq.(|4]), at the minimum; 
with standard values for the constituent quark mass, M = 350 MeV (Mp V /M 2 = 2.25) 
and M = 420 MeV {M PV /M 2 = 1.90), we find, respectively, M N = 1140 MeV and 
M N = 1040 MeV. 

The regularization of the energy, Eq.@, is different from the one used in the calcula- 
tions of refs.|| |J, where in addition to the Dirac continuum also the contribution of the 
bound-state level of the Hamiltonian with the regulator mass was subtracted. The latter 
regularization does not lead to stable solitons. The two regularizations differ by a finite 
contribution of order M 2 /Mp v . From the point of view of computing parton distributions 
both regularizations belong to the class of physically acceptable ones, see below. 

Isosinglet unpolarized quark distribution function. The theoretical framework for the 
calculation of quark distribution functions with the effective chiral theory has been given 
in refs. [||, Q . The isosinglet unpolarized distribution appears in the leading order of the 
l/7V c -expansion and can be represented as a sum of contributions of quark single-particle 
levels in the background pion field: 

u(x) + d(x) = N c f lev (x) + N c J2 [fn(x)-fi°\x)}, (11) 

neg. cont. 

Ux) = M JV |^$t(k)(l + 7 V)^ 3 + K-xM^)$ n (k). (12) 

This expression has been derived from the QCD representation of the parton distribution 
as nucleon matrix element of a light-ray (non-local) operator 0; it is equivalent to the 
definition of the parton distribution as the number of particles carrying a fraction x of the 
nucleon momentum in the infinite-momentum frame ||. The function Eq.(|TT|) describes 
the quark distribution for positive values of x, and minus the antiquark distribution for 
negative x. In Eq.flTTp f^ denotes the matrix element Eq.(|T^) between levels of the 
vacuum Hamiltonian (U = 1), with f^(x) = for x > 0, E n < on kinematical grounds. 
When computing the isosinglet distribution function in the effective theory we must 
consider separately the even and odd part of the function defined by Eg. (|Tl~D , since they 
exhibit different behavior with respect to the UV cutoff. The valence quark distribution 



is given by the even part of Eq. (|Tl"D 



u(x) + d(x) - u(x) - d(x) = N c [f lev (x) + f\ cv (-x)] 

+ N c Y, lfn(x) + f n (-x)-f^(x)-f^(-x)].(l3) 



neg. cont. 



Here the continuum contribution is finite and does not require an UV cutoff. This distri- 
bution is normalized by the baryon number sum rule, 

-l 



dx \u(x) + d(x) - u(x) - d(x)] = N c , (14) 

/o 

which holds "level by level" in the chiral soliton model, as can be seen by integrating 

Eq.(|l2D over x, upon which the 7°7 3 -term drops out. 

The total distribution of quarks plus antiquarks is given by the odd part of Eg. (|TTD . 

The continuum contribution to this function is logarithmically divergent for each value 

of x and requires regularization. (Again, the divergence is proportional to M 2 JJ.) We 

regularize it by a Pauli-Villars subtraction, analogous to the energy, Eq.@: 

u(x) + d(x) + u(x) + d(x) = N c [fi ev (x) - f\ cv (-x)] 

+ N c £ lfn(x) - U-X) - fi°\ X ) + fi°\~x)] 
neg. cont. 

N cj!t- E UpvAx) - IpvA-x) - I ( pUx) + ffU-x)]. (15) 

1V1 PV neg. cont. 

Note that, again, this regularization differs from the one used in refs.|3|, || in that only 
the Dirac continuum is subtracted here but not the level contribution. We emphasize 
that such regularization preserves the general properties of parton distributions just as 
the one used in refs. [§], |J ; in particular, it satisfies all criteria for a "good" regularization 
established in ref.Q], such as absence of anomalies, uniform logarithmic cutoff dependence, 
etc. 

The total distribution of quarks plus antiquarks regularized according to Eq. ([OJ) sat- 
isfies the momentum sum rule, if the classical pion field is a stationary point of the 
regularized energy, Eq.(|)[]: 

/ dx x[u(x) + d(x) + u(x) + d(x)} = 1. (16) 

Jo 

A general proof of the momentum sum rule has been given in ref . || , which can straight- 
forwardly be carried over to the Pauli-Villars regularized energy, Eq.(^j), and distribution 
function, Eq. (|T5|) . Note that the momentum sum rule differs from the baryon number 
sum rule in that it is a statement only about the sum over all occupied levels but does 
not hold "level by level" . 

Results. We compute the distribution functions Eqs. (p~3|, |i~5l) for the self-consistent 



soliton profile using the numerical method developed in ref.[|J. For the numerical cal- 
culations it is convenient to make use of the fact that with Pauli-Villars regularization 



1 We note that for Eq.([16|) to be satisfied it is also necessary that M/v in Eq.([IJ) be the minimum value 
of the static energy, i.e., the nucleon mass in the leading order of the l/iV c -expansion quoted above. 



one may equivalently compute the distribution function at negative x by summing over 
non-occupied states, thus avoiding the need of vacuum subtraction: 

M-x)+ E /»(-*) - E /«(*), (17) 

neg. cont. pos. cont. 

and similarly for the sum with M — > Mpy; see ref.[(|] for details. 

The distribution functions receive contributions from the bound-state level as well 
as from the polarized Dirac continuum of quarks. To understand the nature of these 
contributions it is instructive to consider for a moment the "universal" function, Eq . fllTD , 
which describes the quark distribution at positive and minus the antiquark distribution 
at negative values of x. Fig.[l] shows the contributions of the bound-state level and the 
Dirac continuum to this function (M = 350 MeV), the latter being regularized by a Pauli- 
Villars subtraction as in Eq . flTBl) . The distribution resulting from the bound-state level 
(dashed line) is positive for all values of x; hence it makes a negative contribution to 
the antiquark distribution. The correct sign of the antiquark distribution is obtained by 
including the contribution of the Dirac continuum (dot-dashed line), which is positive for 
x > and negative for x < 0, in agreement with the fact that the logarithmic divergence 
of the distribution function is an odd function x. We note that the so-called valence quark 
approximation advocated in ref.||10||, when applied to parton distribution functions, would 
violate the positivity of the antiquark distributions. 

The total distribution of quarks plus antiquarks, u(x) +d(x) +u(x) + d(x), is obtained 
by taking the odd part in x of the function shown in Fig.|l|, or, equivalently, performing 
the sum Eq. (p~5|) . The resulting distribution is shown in Fig.^| (solid line). It receives 
a sizable contribution from the Dirac continuum (dot-dashed line). We emphasize that 
the calculated distribution explicitly satisfies the momentum sum rule, thanks to the 
self-consistency of the classical pion field. 

The valence quark distribution, u(x) + d(x) — u{x) — d(x), is shown in Fig|| This 
distribution one obtains by taking even part of the function shown in Fig.|l| and removing 
the Pauli-Villars regularization, MpyjM — > oo, or, equivalently, by evaluating Eq. (|l3|) . 
In this distribution the logarithmic divergence of the Dirac continuum cancels; the finite 
continuum contribution obtained in the limit MpyjM — > oo is shown in Fig.§ (dot-dashed 
line); it is numerically small. The calculated distribution satisfies the baryon number sum 
rule, Eq.([TJ|). Also shown in Figs.0 and |3] (triangles) are the distributions calculated in 
ref . pi with Pauli-Villars subtraction of the level contribution and a fixed soliton profile. 

The valence quark distribution receives a non-zero contribution from the Dirac con- 
tinuum (Fig.|3|, dot-dashed line). The contribution of the Dirac continuum to the baryon 
number, Eq. flTj]) , however, is exactly zero. This is in fact necessary for the nucleon 
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to have baryon number unity, since the baryon number of the bound-state level is al- 
ready unity. A Pauli-Villars regularization of the Dirac continuum contribution as in 
u(x) + d(x) + u(x) + d(x), Eq. (|I5"D , would not violate the baryon number sum rule — the 
contributions of the true and the Pauli-Villars continuum to the baryon number would 
be individually zero -- and thus constitute an admissible, if "unnecessary", finite regu- 
larization. Since aside from preserving sum rules and other general requirements the UV 
regularization of the effective theory is fraught with some uncertainty, it is instructive to 
explore the consequences of such finite regularization on the distribution functions. For 
the valence distribution, u(x) + d(x) — u(x) — d(x), the influence of a regularization of 
the Dirac continuum is negligible, since this distribution is dominated by the bound-state 
level contribution, see Fig.^. The antiquark distribution, however, which is defined as the 
difference of Eq.fllBT) and Eq . (|i~3"D , receives a large contribution from the Dirac continuum 
and is potentially more sensitive to finite regularizations of u(x) + d(x) — u(x) — d(x). 
In Fig.[| we show the antiquark distribution obtained as the difference between the same 
distribution of quarks plus antiquarks (Fig.0, M = 350 MeV) and a valence distributions 
computed with i) unregularized continuum contribution (solid line) and ii) continuum 
Pauli-Villars regularized with Mp V /M 2 = 2.52 as in the total distribution (short-dashed 
line). One sees that the effect of the finite regularization on the antiquark distribution 
is not dramatic. In particular, the antiquark distribution is positive in both cases. Of 
course, it is guaranteed to be positive in the limit M PV /M — > oo, because of the log- 
arithmic divergence of u(x) + d(x) + u(x) + d(x). The fact that it is positive also for 
finite Mpy/M shows that the finite, regularization-dependent terms do not upset this 
asymptotic statement for the particular regularization used here, as it should be. 

Figs.0 and [^ also show the distributions obtained for the soliton with M = 420 MeV 
(dotted lines). We see that the distributions calculated with the self-consistent soliton 
profile are rather insensitive to the value of the constituent quark mass. This parallels 
the behavior of most hadronic observables in this model Jll| . 



The numerical method employed here |4| does not allow for an accurate determination 
of the limiting behavior of the parton distributions for x —>■ 0. The limiting behavior of 
the Dirac continuum contribution in Fig.|l| has been obtained by extrapolation, see ref . [§j 
for details; for the distributions multiplied by x, Figs]^, ^ and SI this limitation of the 
present method is irrelevant. We note also that the "small-x" behavior of the parton 
distributions in the effective chiral theory, i.e., the behavior in the parametrically small 
region \x\ < (M/A) 2 /iV c , depends on the details of the UV regularization, as has recently 



been emphasized in connection with off-forward parton distributions |12 



The parton distributions studied here pertain to the large- N c limit, where the nucleon 



is heavy, M^ oc N c . The calculated distributions therefore do not go to zero at x — 1, 
rather, they are exponentially small at large x, as discussed in ref.[[§. 

Summary and Conclusions. We have shown, by a specific example, that the general 
requirements on the UV regularization in the calculation of parton distributions in the 
effective chiral theory can be reconciled with the requirements of a stable soliton solution, 
which depends on 0(M 2 /A 2 )-terms in the energy functional. The calculated quark and 
antiquark distributions explicitly satisfy the baryon number and momentum sum rules. 

The distributions obtained with the self-consistent soliton profile largely support the 
results of the calculations of refs. H §] using a fixed soliton profile. Except for the momen- 
tum sum rule the differences are not significant, given the limited theoretical accuracy of 
this approach. We note that the regularization used in the present work is only an exam- 
ple of a class of possible regularizations leading to a stable soliton; for instance, one may 
perform multiple Pauli-Villars subtraction of the Dirac continuum with different masses, 
etc. However, we see at present no need to explore such possibilities. 

The basic expressions for the parton distributions in refs.jS], [| have been derived 
in the leading order of M/A, in which the "constituent" quarks of the effective theory 
can be regarded as pointlike. This is why the momentum sum rule is saturated by the 
quark and antiquark distributions computed here. If one knew the precise form of the 
effective theory at order M 2 /A 2 , one could "resolve" the structure of the constituent 
quark and recover the distribution of quarks and gluons inside the constituent quark. 
Such a picture is in fact suggested by the instanton vacuum, where the parameter M 2 /A 2 



is related to the small packing fraction of the instanton medium [Jl^, [14]; its concrete 
realization is the subject of current investigations. It is not unlikely that the "resolution" 
of the constituent quark structure could be represented in the form of a convolution 
in Bjorken-x of the constituent quark distributions with parton distributions inside the 
constituent quark. This general idea is in fact rather old |15|| ; it has recently been taken 



up again fll6| . Since the "constituent" quark and antiquark distributions computed here 



are correctly normalized to the total nucleon momentum, they can serve as a starting 
point for investigations assuming a structure inside the constituent quark. 
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Figure 1: The isosinglet distribution function, Eq.flTID, corresponding to u(x) + d(x) 
for x > and — u{— x) — d(—x) for x < 0, for the self-consistent soliton with M = 
350 MeV. Dashed line: contribution of the bound-state level, f\ ev (x), giving a negative 
contribution to the antiquark distribution. Dot-dashed line: contribution of the Pauli- 
Villars regularized Dirac continuum. Solid line: total result (bound-state level plus Dirac 
continuum). The total antiquark distribution is positive. 
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Figure 2: The isosinglet distribution of quarks plus antiquarks, \x[u(x) + d(x) + u(x) + 
d(x)]. Shown are the total result (solid line) and the contribution of the Dirac continuum 
(dot-dashed line) for the self-consistent soliton with M = 350 MeV. Also shown is the 
total result for the soliton with M = 420 MeV (dotted line). Triangles: the distribution 
calculated in ref.pj with an analytic soliton profile and regularized level contribution. 
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Figure 3: The isosinglet valence quark distribution, ^x[u(x) + d(x) — u(x) — d(x)]. Shown 
are the total result (solid line), and the contribution of the unregularized Dirac continuum 
(dot-dashed line) for the self-consistent soliton with M = 350 MeV. This distribution is 
obtained from the symmetric part of the function in Fig.|l] after removing the Pauli-Villars 
cutoff, Mpv/M — > oo. Also shown is the total result for the soliton with M = 420 MeV 
(dotted line). Triangles: the distribution calculated in ref.[|] with an analytic soliton 
profile. 
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Figure 4: The isosinglet antiquark distribution, x[u(x) + d(x)], for the self-consistent 
soliton with M = 350 MeV, obtained as difference between u(x) + d(x) + u(x) + d(x), 
Fig.|2], and u(x) + d(x) — u(x) — d(x), Fig.[| with different regularizations of the continuum 
contribution to u(x) + d(x) — u(x) — d(x). Solid line: unregularized. Short-dashed line: 
regularized with M 2 PV /M 2 = 2.25. 



13 



